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ABSTRACT 

Aims. The aim of the project is to improve our knowledge on the low-mass and low-metallicity population to investigate the influence 
of metallicity on the stellar (and substellar) mass function. 

Methods. We present the results of a photometric and proper motion search aimed at unearthing ultracool subdwarfs in large-scale 
surveys. We employed and combined the Fifth Data Release (DR5) of the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area 
Survey (LAS) and the Sloan Digital Sky Survey (SDSS) Data Release 7 complemented with ancillary data from the Two Micron 
All-Sky Survey (2MASS), the DEep Near-Infrared Survey (DENIS) and the SuperCOSMOS Sky Surveys (SSS). 
Results. The SDSS DR7 vs UKIDSS LAS DR5 search returned a total of 32 ultracool subdwarf candidates, only two being recognised 
as a subdwarf in the literature. Twenty-seven candidates, including the two known ones, were foUowed-up spectroscopically in the 
optical between 600 and 1000 nm thus covering strong spectral features indicative of low metallicity (e.g., CaH) , 21 with the Very 
Large Telescope, one with the Nordic Optical Telescope, and five were extracted from the Sloan spectroscopic database to assess (or 
refute) their low-metal content. We confirmed 20 candidates as subdwarfs, extreme subdwarfs or ultra-subdwarfs with spectral types 
later than M5; this represents a success rate of >60%. Among those 20 new subdwarfs, we identified two early-L subdwarfs very likely 
located within 100 pc that we propose as templates for future searches because they are the first examples of their subclass. Another 
seven sources are solar-metallicity M dwarfs with spectral types between M4 and M7 without Ha emission, suggesting that they 
are old M dwarfs. The remaining five candidates do not have spectroscopic follow-up yet; only one remains as a bona-fide ultracool 
subdwarf after revision of their proper motions. We assigned spectral types based on the current classification schemes and, when 
possible, we measured their radial velocities. Using the limited number of subdwarfs with trigonometric parallaxes, we estimated 
distances ranging from ~95 to ~600 pc for the new subdwarfs. We provide mid-infrared photometry extracted from the WISE satellite 
databases for two subdwarfs and discuss their colours. Finally, we estimate a lower limit of the surface density of ultracool subdwarfs 
of the order of 5000-5700 times lower than that of solar-metallicity late-M dwarfs. 

Key words. Stars: subdwarfs — techniques: photometric — techniques: spectroscopic — Infrared: Stars — surveys — Virtual 
Observatory tools 



1 . Introduction adopted classification for M subdwarfs (sd M) and extreme sub - 
H^,,, ^ , , n ■ ,.TT, r ,. , dwarfs (esdM) has recently been revised bv lLepine et al.l(l2007h . 
_ _ Cool subdwarfs ai-e metal-deficient population II dwarfs which ^ ^j^^^ subdwarfs, the ultra-subdwarfs (usdM), has been 
appeal- less luminous than their solai-metallicity counterparts ^^^^^ ^^e sdM and esdM classes originally defined bv lGi^ 
dueto the dearth of metals in their atmospheres ( |Baraffe et al. | ^^^^^^ ^^^^^ ^ parameter, -TT.o/CaH, which 
imji,. They tend to exhibit halo or thick disk kinematics in- „tifle, ^^e weakening of the strength of the TiO band (in 
cluding^noticeable proper motion and large heliocentric yeloc- ^j^^ -^^^^ ^ ^^^^^-^^ metalUcity. An alternative classi- 
ities j Gizis , 1997) . They are very ol d and represent usefu trac- ^^^^.^^ ^^^^^ temperature, gravity, and metallicity has been 
ers of the Galactic chemic al history ( |Burgasseret al. || 2003| ). The p.^posed by Jao et al. ( 2008). The range of metallicity for sub- 
Send ojfprint requests to: N. Lodieu dwarfs, exti'eme subdwaifs, and ultra-subdwarfs span appro xi- 
* Based on observations made with ESO Telescopes at the La Silla mately -0.5 and -1.0, -1.0 and -1^ and below -1.5, respec- 
Paranal Observatory under programme ID 084.C-0928A lively (|Gizis||1997|; |Woolf et al.||2009|). M-type subdwarfs have 
** Based on observations made with the Nordic Optical Telescope, op- typically eff^ective temperatures below -3500-4000 K (depend- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, mg on the metalhcity Baraff^e et al. 1997; Woolf et al. 2009) and 
Norway, and Sweden, in the Spanish Observatorio del Roque de los should display high gravity (log g~5. 5) although sorne variations 
Muchachos of the Instituto de Astroflsica de Canarias. is seen among low-metallicity spectra jJao et al.ll2008l) . 
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Subdwarfs were generally identified from optical (B j, R, and 
/) proper motion catalogues on photographic plates at different 
epochs (Luyten 1979, 1980; Scholz et al. 2000; Lepine et afl 
I2003ht fLodieu et alJ |2005). Several surveys have been con- 
ducted to searc h for subdwar fs over a large temperature range, 
including hot (RvanI 19891) . intermediate (Yong & Lambert 
l2003t iDigbv et al.1 l2003h . and cool components ( Gizisi 11997^ 
A growing number of M subdwarfs have been announced over 
the past years, raising the number of metal-defi cient d warfs 

(iGizisI 
etail 



ten 



with spectral types later than M7 to about 
Gizis & Reid 1997; Schweitzer et al. 1999; Lepine 



1997 



l.ll2003a: 



Scholz et alj2 004b; Gushin g et al.l200 9). This number increased 



significantly after the discovery of 23 late-type subdwarf s in the 
Sloan Digital Sky Survey (SDSS; Adelman-McCarthv & et all 
|2009|) spectroscopic database (Le pine & Scholzi i2008 l) and 
15 others in the multi-epoch database of the 2MASS survey 
dKirkpatrick et al.l 1201 0). Moreover, the hydrogen-burning 
limit has been crossed with the discov ery of the first substella r 
subdwai-f in the 2MASS database (iBurgasser et al.l l2003h . 
This discovery was quickly followed by the announcement 
of another L subdwarf (Burgasser 2004) and, more recently, 
by seven new ones (jSivarani et al. 2009; Gushing et al. 200! 



Lodieu et al] |20 1 Ol; "^ Kh -kpatrick et al. .2010; ,Schmidt et al 



2010al: Bowler et al ] l2010r This number of ultracool subdwarfs 



remains however very small and is at odds with the numerous 
L and T dwarfs reported in the solar neighbourhood within the 
framework of large-scale optical and infrared s urveys, including 
the Two Micron All Sky Surve y (2MASS; e.g. lKirkpatrick et all 



120001; IBurgasser etalT 120021), DEep Near Infrared Survey 
(DENIS; e.g. Delfosse et al. 1997", 49991 'Martin et al." 1999), 
SDSS (e.g. Fanetal. 2000; Leggettetal. 2000; Geballe et ak 
200l. the UKIRT I nfrared Deep Sky Su rv ey (UKIDSS; e.g 



Lodieu et alj l2007t iPinfield et alj 120081: iBumingham et all 
the Ganada-France-Hawaii Brown Dwarf s urvey 



dDelorme et al.' 2008"; Revle et al. 2010; Albert et alj|2011|), and 
WISE CKirkpatrick et aL,201 k.Gushing et al..,201 li) . 



The aim of our study is to identify a complete census 
of metal-poor dwarfs to bridge the gap between the coolest 
M subdwarfs and the recent L subdwarfs identified in vari- 
ous surveys. In this paper we report the outcome of a photo- 
metric and proper motion se arches using the UKI RT Infrared 
Deep Sky Survey (UKIDSS: ITawrence et al]l2007h Large Area 
Survey (LAS) Data Release 5 (DR5) an d the Sloan Digital Sky 
Survey (SDSS) Data Release 7 (DR7; lAbazaiian et al.l l2009l) . 
This catalogu e search was complemented by ancil lary data 
from 2MASS (iGutri et al.ll2003t ISkrutskie et al.ll2006l). DENIS 
(IDENIS Consortiumll2005h . and SuperGOSMOS ("Hamb lv et all 
[2001c b a|). In Sect. |2] we describe the photometric and proper 
motion criteria designed to identify ultracool subdwarfs in pub- 
lic databases using Virtual Observatory tools. In Sect.^we detail 
the spectroscopic follow-up conducted at optical wavelengths 
with the visual and near UV FOcal Reducer and low dispersion 
Spectrograph (FORS2; Appenzeller et al. 1998) mounted on the 
Very Large Telescope (VLT), the ALFOSG spectrograph on the 
Nordic Optical Telescope (NOT) and complemented by opti- 
cal spectra downloaded from the SDSS spectroscopic database. 
In Sect. |5] we present the results of the spectroscopic analysis 
and determine the main properties of the newly identified ultra- 
cool subdwarfs including colours, spectral types, radial veloci- 
ties, and distances. Finally, we summarise our work in Sect. |6] 
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Fig. 1. {J - Ks,i - J) colour-colour diagram showing the 
positions of our subdwarf candidates identified in the 
SDSS/UKIDSS search (filled circles) with their respective er- 
ror bars. Overplotted as filled red triangl e s are two known ul- 
tracool subdwarfs from iLepine & Scholzl (|2008|) with J and K 
photometry from UKIDSS LAS DR5 (error bars are smaller 
than the symbols) and the solar-metallicity M/L dwarf sequence 
(line with crosse s; spectral t ypes are labell ed; OsMO, lOsLO) 
from IWest etaP (I200&) and Ischmidt et al.1 (HoiOb). The pho- 
tometry shown on t his diagram in the MKO photometric system 
dHewett et al.ll2006h . 



2. Sample selection 

In this section we describe the selection procedure developed to 
unveil new subdwarfs in the cross-correlation of the SDSS DR7 
and UKIDSS LAS DR5 catalogues. 

The idea of a generic photometric search for subdwarfs was 
triggered by the (/ - J,J - K) colour-colour diagram is pre- 
sented in Fig. 4 of I Scholz et al . (2004b). A similar diagram 
is shown in Fig. [Tl Kn own ultracool subdwarfs identified by 
^Lepine & Scholzl (l2008h in the SDSS spectroscopic d atabase are 
plotted al ong with the sequ ence of field M/L dwarfs (IWest et al.l 
2008; Schmidt et al. 2010b). Subdwarfs tend to be bluer than so- 
lar metallicity M/L dwarfs but occupy the same region in the i-J 
domain as their solar metallicity counterparts, making their iden- 
tification difficult in a pure optical-to-infrared search. However, 
the onset of collision-induced molecular hydrogen H2 opacity 
at near-infrared wavelengths generates bluer J - H and J - Kg 
colours {J - Ks < 0.7) than observed for solar metallicity M and 
L dwarfs. 

To optimise our photometric selection, we have employed 
the reduced proper motion (Hr=r -H 5x log(/i) -1- 5) as a proxy 
for metallicity. This parameter is useful to separate solar- 
metal licity s tars from subdwarfs, and white dwarfs ( Tone si 



197 ilEvansiri9 92: Salim & Gould 20(^ lLepine & Sharalf2005r 
Burgasse r et al.l2007i:iLodieu et al. 200y). A reduced proper mo- 
tion diagram is s hown in Fig. |2] with known subdwarfs from 
Lep ine & Scholzl (|2008) marked as open symbols and our new 
candidates confirmed as subdwarfs shown as filled symbols. 
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Table 1. Running number, names according to the lAU nomenclature of UKIDSS, coordinates (in J2000), SDSS optical (ugriz) magnitudes with 
epoch, UKIDSS near-infrared (YJHK) magnitudes with epochs, total proper motion (yu in arcsec/yr), and the reduced proper motion (H,) of 
ultracool subdwarf candidates extracted from the cross-match between the UKIDSS LAS DR5 and SDSS DR7 



ID 



Name 



R.A. 
hh:mm:ss.ss 



Dec 



SDSSu 
mag 



SDSSg 
mag 



SDSSr 
mag 



SDSS/ 
mag 



SDSSz 
mag 



Epoch 



mag 



mag 



H 

mag 



K 
mag 



Epoch 



/J 

"/yr 



mag 



1 ULAS J004539.97+135032.7 

2 ULAS J012830.89+134507.4 

3 ULAS JO 1 3346.25+ 1 32822.4 

4 ULAS J015034.33+142002.4 

5 ULAS J02 533.75+123824.1 

6 ULAS J033350.84+001406.1 

7 ULAS J083525.4I-000940.5 

8 ULAS J084153. 89+020615. 1 

9 ULAS J084358.50+060038.6 

10 ULAS J085500. 12+000204.1 

11 ULAS J093502. 10+8 551.8 

12 ULAS J095133. 35+070208.9 

13 ULAS J100126.29-013426.6 

14 ULAS J100743.96-022830.0 

15 ULAS J101613. 89+11311.4 

16 ULAS Jl 15826.62+044746.8 

17 ULAS J120214.62+073113.8 

18 ULAS J1215 8.37+040200.5 

19 ULAS J122145.28+080404.4 

20 ULAS J123659.43-002158.2 

21 ULAS J124234.62+143306.2 

22 ULAS J 1 24425 . 90+ 1 0244 1 . 9 

23 ULAS J124621.90+044309.9 

24 ULAS J125635.91-001944.9 

25 ULAS J131705. 66+091016.9 

26 ULAS J141806.71+000035.5 

27 ULAS J145441.42+123556.7 

28 ULAS J15121 1.64+064251.3 

29 ULAS J 15433 1 .93+24537.8 

30 ULAS J154450.60+051613.6 

31 ULAS J154855. 74+080508.2 

32 ULAS J233359.39+04935.2 



00:45:39.97 
01:28:30.89 
01:33:46.25 
01:50:34.33 
02:05:33.75 
03:33:50.84 
08:35:25.41 
08:41:53.89 
08:43:58.50 
08:55:00.12 
09:35:02.10 
09:51:33.35 
10:01:26.29 
10:07:43.96 
10:16:13.89 
11:58:26.62 
12:02:14.62 
12:15:08.37 
12:21:45.28 
12:36:59.43 
12:42:34.62 
12:44:25.90 
12:46:21.90 
12:56:35.91 
13:17:05.66 
14:18:06.71 
14:54:41.42 
15:12:11.64 
15:43:31.93 
15:44:50.60 
15:48:55.74 
23:33:59.39 



+ 13:50:32.7 
+ 13:45:07.4 
+ 13:28:22.4 
+ 14:20:02.4 
+ 12:38:24.1 
+00:14:06.1 
-00:09:40.5 
+02:06:15.1 
+06:00:38.6 
+00:02:04.1 
+08:05:51.8 
+07:02:08.9 
-01:34:26.6 
-02:28:30.0 
+01:13:11.4 
+04:47:46.8 
+07:31:13.8 
+04:02:00.5 
+08:04:04.4 
-00:21:58.2 
+ 14:33:06.2 
+ 10:24:41.9 
+04:43:09.9 
-00:19:44.9 
+09:10:16.9 
+00:00:35.5 
+ 12:35:56.7 
+06:42:51.3 
+02:45:37.8 
+05:16:13.6 
+08:05:08.2 
+00:49:35.2 



24.559±1.069 
24.930± 1.485 
23.464±0.835 
23.456±0.617 
23.281±0.597 
23.392±0.753 
24.659±0.781 
23.580±0.904 
22.352±0.295 
24.020±0.725 
23.088±0.548 
25.254±0.829 
25.659±0.956 
25.722±0.798 
24.805±1.451 
24. 243 ±0.864 
24.028±0.953 
24.584±0.971 
26.002±0.968 
24.089±0.876 
25.128±0.972 
25.453±0.726 
23.343±0.639 
24.973±0.772 
25.781±0.946 
25.359±0.824 
22.972±0.350 
25.525±1.651 
25.682±0.934 
23.830±1.103 
24.355±1.431 
24.651 ±0.928 



25.657+0.771 
26.014+0.640 
20.822+0.035 
22.912+0.196 
21.944+0.104 
23.818+0.374 
24.601+0.532 
26.662+0.390 
19.711+0.015 
24.205+0.400 
24.742+0.709 
24.666+0.677 
24.326+0.700 
25.124+0.871 
24.590+0.988 
23.832+0.338 
22.681+0.130 
23.636+0.381 
26.215+0.756 
21.392+0.050 
21.662+0.059 
25.472+0.566 
21.912+0.095 
21.414+0.045 
22.569+0.147 
24.845+1.117 
21.191+0.035 
24.856+1.050 
25.591+0.715 
24.671+0.770 
25.155+0.720 
23.284+0.200 



21.500±0.075 
22.595±0.294 
18.955±0.012 
21.087±0.058 
19.767±0.021 
21.731±0.102 
22.497±0.151 
23.681±0.697 
17.815±0.006 
22.320±0.128 
22.374±0.185 
22.335±0.178 
22.214±0.170 
23.240±0.440 
22.513±0.260 
21.961±0.093 
20.601±0.034 
21.517±0.063 
22.897±0.430 
19.563±0.015 
19.794±0.020 
22.499±0.157 
20.083±0.024 
19.484±0.014 
20.505±0.038 
21.872±0.102 
19.274±0.012 
22.424±0.284 
23.278±0.383 
22.769±0.192 
22.719±0.214 
21.476±0.066 



20.254±0.040 
21.154±0.121 
17.822±0.008 
19.645±0.026 
18.113±0.010 
19.193±0.017 
20.680±0.047 
21.821±0.253 
16.717±0.005 
20.745±0.048 
21.318±0.121 
20.994±0.092 
21.051±0.107 
21.528±0.142 
21.036±0.113 
19.426±0.017 
19.444±0.018 
20.016±0.028 
21.635±0.227 
18.404±0.010 
18.740±0.014 
19.483±0.020 
19.000±0.016 
18.231±0.009 
19.313±0.017 
20.286±0.041 
18.107±0.009 
21.189±0.118 
21.929±0.187 
21.449±0.114 
21.178±0.070 
20.152±0.036 



19.524+0.079 
20.447+0.256 
17.181+0.013 
18.895+0.040 
17.287+0.018 
17.811+0.021 
19.851+0.088 
20.290+0.222 
16.083+0.008 
19.929+0.107 
20.494+0.197 
20.001+0.131 
20.444+0.206 
20.461+0.196 
20.514+0.261 
18.155+0.022 
18.720+0.038 
19.116+0.047 
20.358+0.235 
17.745+0.018 
18.121+0.028 
18.019+0.019 
18.415+0.041 
17.518+0.015 
18.561+0.032 
19.522+0.085 
17.449+0.013 
20.351+0.208 
21.047+0.366 
20.978+0.407 
20.786+0.241 
19.427+0.063 



1999.78 
1999.78 
1999.78 
1999.78 
2001.89 
1999.78 
1999.22 
2000.92 
2002.19 
2000.17 
2002.19 
2002.93 
2000.17 
2000.17 
2000.98 
2001.14 
2003.25 
2001.14 
2003.25 
1999.22 
2003.08 
2002.19 
2001.29 
1999.22 
2003.32 
1999.22 
2003.47 
2003.32 
2000.36 
2001.46 
2003.25 
2001.79 



18.755±0.055 
19.693±0.172 
16.351±0.009 
18.018±0.037 
16.456±0.008 
16.814±0.006 
19.098±0.057 
19.638±0.094 
15.268±0.003 
19.293±0.072 
19.767±0.108 
19.384±0.063 
19.517±0.088 
19.891±0.098 
19.573±0.111 
17.121±0.014 
18.059±0.046 
18.302±0.039 
19.973±0.087 
16.900±0.016 
17.334±0.017 
16.982±0.008 
17.664±0.018 
16.694±0.009 
17.843±0.025 
18.595±0.055 
16.664±0.007 
19.501±0.078 
19.536±0.092 
19.685±0.117 
18.783±0.045 
18.541±0.058 



18.043+0.059 
19.318+0.155 
15.855+0.010 
17.424+0.032 
15.872+0.009 
16.106+0.010 
18.259+0.077 
18.995+0.097 
14.732+0.004 
18.610+0.072 
19.182+0.108 
18.891+0.083 
18.779+0.082 
19.060+0.130 
18.855+0.146 
16.486+0.016 
17.539+0.029 
17.698+0.039 
19.141+0.115 
16.406+0.012 
16.821+0.019 
16.259+0.010 
17.141+0.020 
16.167+0.009 
17.311+0.027 
18.079+0.031 
16.121+0.010 
18.830+0.080 
18.742+0.099 
19.232+0.108 
18.214+0.050 
18.018+0.036 



17.630±0.088 
18.585+0.145 
15.467±0.015 
16.946+0.042 
15.709+0.012 
15.765+0.009 
18.186+0.099 
18.496+0.114 
14.307+0.004 
18.366+0.094 
18.685+0.109 
18.311+0.088 
18.469+0.117 
18.534+0.106 
18.447+0.115 
16.124+0.014 
17.041+0.033 
17.242+0.046 
18.417+0.094 
16.040±0.015 
16.485±0.023 
15.999±0.011 
16.718±0.024 
15.729±0.011 
16.940±0.034 
17.723±0.090 
15.695±0.012 
18.503±0.109 
18.160±0.124 
18.604±0.127 
17.745+0.083 
17.628±0.070 



17.363±0.100 
18.719±0.301 
15.250±0.017 
16.794±0.069 
15.590±0.018 
15.504±0.019 
17.917±0.138 
18.304±0.194 
14.052±0.006 
18.344±0.211 
18.527±0.194 
18.374±0.175 
18.198±0.199 
18.398±0.227 
18.451±0.231 
15.836±0.020 
16.902±0.051 
17.154±0.076 
18.488±0.217 
15.775±0.023 
16.283±0.036 
15.773±0.019 
16.447±0.035 
15.518±0.017 
16.738±0.049 
17.801±0.183 
15.471±0.013 
18.181±0.160 
18.049±0.194 
18.693±0.214 
17.633±0.127 
17.368±0.090 



2007.71 
2007.64 
2007.63 
2007.64 
2006.89 
2007.63 
2006.93 
2006.99 
2006.96 
2006.93 
2007.06 
2008.13 
2005.99 
2005.97 
2008.13 
2008.13 
2007.29 
2007.26 
2007.06 
2005.44 
2007.26 
2007.12 
2007.05 
2005.45 
2006.47 
2005.45 
2007.26 
2006.43 
2006.57 
2005.43 
2006.44 
2007.61 



0.14 
0.13 
0.28 
0.25 
0.27 
0.78 
0.17 
0.30 
0.50 
0.16 
0.61 
1.01 
0.80 
0.27 
0.16 
0.19 
0.34 
0.25 
0.30 
0.29 
0.29 
0.67 
0.18 
0.36 
0.32 
0.19 
0.31 
0.33 
0.73 
1.21 
1.13 
0.20 



22.177 
23.240 
21.211° 
23.106 
21.925 
26.201 
23.592 
26.044 
21.307* 
23.364 
26.314 
27.353 
26.739 
25.421 
23.556 
23.346 
23.284'" 
23.545 
25.306 
21.877'' 
22.133 
26.638 
21.334 
22.236'' 
23.030 
23.248 
21.742^' 
24.992 
27.602 
28.182 
27.983 
23.005* 



Notes: " LP 46 8-277 jLepine & SharJ2005HWest et al.l2008ft 

* LHS 2045 (ID=9) I'Lepine & Shara"2005!;'West et al.'2008') 
" SDSS J 1202 14.6 1+0731 13.7 (ID=17) classified as M4 by West etani200a> 
'' 2MASS J12365942-0021578 (1D=2) classified as M3 by West eV^UM^ 
' 2MASS J12563590-00194507 (ID=24) classified as M3 by We st et alj 1200^) 

I 2MASS J14544145+1235576 (lD=27: .LeDiiie & Shara.2005.) 

s Located at -42 arcsec from 2MASS J23335840+00501 19 classified as a LO dwarf bv lztiang et al.l OoToh 
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Our search has been conducted taking advantage of 
Virtual Obsei-vatory (VOfl tools, namely STILTSEI and Aladir£| 
(Bonnarel et al. 2000). The detailed photometric and proper mo- 
tion criteria used for the SDSS DR7 vs UKIDSS LAS DR5 
search are given below. The resulting total number of candidates 
returned by this query is 33 (Table [1} but one was rejected af- 
ter looking at the images. Those candidates are shown as filled 
symbols in Figs.[T]and|2] The finding charts created with Aladin 
are displayed in Figs. [T] and [2] of the Appendix. 

- Only SDSS point sources (class=6) were considered 

- UKIDSS point sources (mergedClass parameter equal to 
-1 or -2) and detections in Y, J, and H 

- Only sources fainter than Y > 10.5, J > 10.5, and H > 10.2 
mag were considered to avoid saturated stars 

- For all SDSS sources, we looked for UKIDSS counterparts 
between 1 and 5 arcsec, implying that we are sensitive to 
dwarfs with proper motions between 0.125 and ~1.0 arcsec 
per year depending on the baseline existing between SDSS 
and UKIDSS (our targets span the 3.1 1-7.93 year baseline) 

- Colour selection were then applied as follows: r - i > 1.0 
mag, g - r > 1.8 mag, r - z > 1.6 mag 

- Only UKIDSS sources with good quality flags were kept i.e. 
ppErrBits parameter in each filter less or equal to 256 

- Infrared criterion of J - K < 0.7 mag was applied 

- Xi and Eta parameters refering to positional matching 
should be between -0.5 and 0.5 

- Only objects satisfying H, > 20.7 mag were kept 

The proper motion of each candidate was computed using 
the UKIDSS LAS and SDSS positions only but images from 
2MASS, DENIS, and SuperCOSMOS were checked by eye for 
additional epochs to get rid off false positives. A revision of the 
proper motions is detailed in Section 1X2] 

Candidates were visually inspected using the scripting ca- 
pabilities of Aladin. Sources from DENIS, 2MASS, SDSS, 
SuperCOSMOS and UKIDSS as well images from UKIDSS, 
2MASS and SDSS were used in the analysis. False candidates 
were rejected due to several reasons, the most likely being the 
mismatch between the SDSS source and the associated UKIDSS 
counterpart. 

We imposed a lower limit on the proper motion on pur- 
pose in order to bias our search towards halo objects (e.g., 
IScholz et alj2000l) and avoid contamination by "normal" dwarfs 
and extragalactic sources. Originally, we took 0.5 arcsec (in- 
stead of 1 arcsec) as the lower limit for the SDSS-UKIDSS 
separation but it produced a large number of false candi- 
dates. Moreover, imposing a detection in the g-band is ham- 
pering the detection of cooler subdwaifs which would be 
faint at short optical wavelengths. We tested our complete- 
ness limit by checking how many of the ~250 known ultra- 



work (iGizisI 19971 


Schweitzer et al. 1999 


LeDine et al.' '20021 


2003bl; Scholz et al. 


I2004blla: Lodieu etal 


2005 


Reid & Gizis 


20051 iLepine & Scholz* '2008': 'Kirkpatrick 


et al. 


2010). Of the 



~250 M- and L-type subdwarfs, only 11 lie in the common 
area between UKIDSS LAS DR5 and SDSS DR7: seven do 
not satisfy either the optical colour constraints because they 
have spectral types earlier than M5 or are missing near-infrared 



' The Virtual Observatory (|http://vyvy w.ivoa.rietl l is an operational re- 
search infrastructure designed to facilitate the access and analysis of the 
information hosted in astronomical archives 

^ http://www.star.bris.ac.uk/~mbt/stilts/ 

^ http://aladin.u-strasbg.fr/ 



photometry in at least one of the three bands (YJH); two 
out of 11 are recovered by our search criteria: LHS2045 
(esdM4.5) and S DSS J020533.75-I-123824.0 (ID=5; esdM8.5; 
iLepine & Schoi3 l2008g). However, two a re not within our sam- 
ple: LHS2096 (e sdM5.5; Hr=20.62 mag ILepine & Sharall200l 
IWest et all l2008l) because of the lower limit we imposed on 
the reduced proper motion, and SDSS J023557.61 -1-010800.5 
(esdM7; Lepine & Scholz 2008) because of its separation be- 
tween the SDSS and UKIDSS coordinates (0.92 arcsec) com- 
pared to our lower limit of 1.0 arcsec. 



3. Proper motion revision 

In this section we discuss the accuracy of the proper motions 
obtained by the Virtual Observatory tools by comparing the po- 
sitions in the SDSS and UKIDSS catalogues (no errors consid- 
ered) and dividing by the epoch difference. 



3.1. Method 

We computed accurate proper motions by measuring the pixel 
coordinates {x,y) of the targets and tens of other sources on the 
Sloan z and UKIDSS LAS J images. We carried out this proce- 
dure under the IRAF environment. Firstly, we downloaded the 
SDSS and UKIDSS LAS images with a size of 6 arcmin aside, 
centered on each target. Secondly, we identify high-quality point 
sources with a signal-to-noise higher than 10 in both images, ob- 
jects selected as reference stars for proper motion measurements. 
We assumed that these sources (about 30 per target) are not mov- 
ing, assumption valid because they are centered around (0,0) in 
Fig.S 

Then, we transformed the pixel coordinates from one epoch 
to the other epoch using second-order polynomial transforma- 
tions for both X and y axes. The dispersion of the transforma- 
tions was typically 0.16 pixel (or 0.032 arcsec). The resulting x 
and y shifts were converted into proper motions by taking into 
account the time baseline of the data and the appropriate pixel 
scale values. Our measurements along their associated error bars 
are given in Table |2] 

3.2. New proper motions 

We find that the revised proper motions agree with the proper 
motion obtained from the VO for 22 of the 32 candidates (~70%) 
(Table |2]i Below we give details on the 10 sources whose new 
proper motion differ from the original ones, leading to a different 
position in the reduced proper motion diagram (Table |2] Fig. |2]l: 

- Sources #30 and #3 1 should be rejected because no object 
is detected on the UKIDSS images suggesting that these are 
false cross-matches. No optical spectra (Section|4]i are avail- 
able for these sources so we reject them from our sample 

- Sources #11 and #12 have smaller proper motion because 
of large offsets between the position of the source on the 
SDSS image and the position reported by the SDSS cata- 
logue. Their revised Hr is smaller than 20.7 mag, implying 
that they would not enter our sample. No optical spectra are 
available for them so we reject them 

- Three sources classified as M dwarfs (#8, 28, and 29) are 
rejected because their revised proper motions are smaller 
than the one derived by the VO, leading to smaller Hr value 
and/or a position in the reduced proper motion diagram sug- 
gesting that they are solar-metallicity dwarfs 
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Fig. 2. Reduced proper motion diagram for new ultracool subdw arfs identified in the S PSS D R7 vs UKIDSS LAS DR5 search. 
The small dots represent all sources in the catalogue published by Lepine & Sharal (LS05;|2P05) with counterpart in the SDSS DR7 
database. Our new subdwarfs, extreme subdwarfs, and ultra- subdwarfs are marked as filled squares, circles, and triangles, respec- 
tively. Diamonds represent our candidates classified as solar- me tallicity M dwarfs. The five sources with out optical spectroscop y 
are marked with a star symbol. Known subdwarfs from the literature are marked as open symbols (LS08: lLepine & Scholz 2008!). 
Left: Diagram for the values given by the Virtual Observatory. Right: Diagram using the revised proper motions computed with the 
method detailed in Section[3] 
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Fig. 3. Revised proper motions of our candidates. Astrometric 
reference stars (small open squares) are concentrated around 
(0,0). Symbols as in Fig.|2] Error bars are smaller than the size 
of the symbols. 

- Three M dwarfs (#14, 15, and 19) have smaller revised 
proper motions but lie in the same region as subdwarfs in 
the reduced proper motion diagram 



We note that sources #7 and #2 have revised proper motions 
identical to those from the VO and both objects lie in the subd- 
warf domain in the reduced proper motion diagram. Source #7 
has no optical spectroscopy, it thus remains a reliable ultracool 
subdwarf candidate. On the contrary, source #2 has optical spec- 
troscopy suggesting it is a solar-metallicity M dwarf contaminat- 
ing our sample. 

To summarise, we conclude that only one of the five sources 
without optical spectrum remains as a bona-fide subdwarf. 
Moreover, three of the seven sources spectroscopically classi- 
fied as M dwarfs (Fig.|6]l would be rejected based on the revised 
proper motion, the remaining ones being photometric contami- 
nants because they lie in the region of the reduced proper motion 
diagram where ultracool subdwarfs are found (right panel in Fig. 

4. Optical spectroscopy 

We emphasise that optical spectroscopy was obtained for can- 
didates identified in the preliminary search i.e. with the proper 
motion derived from the VO and not from the revised proper 
motions discussed in Section |3] 

4.1. NOT spectroscopy 

We carried out low-resolution (R~450) optical (500-1025 nm) 
spectroscopy with the ALFOSC spectrograph on the 2.5-m 
Nordic Optical Telescope (NOT) at the Obsei-vatory of the 
Roque de Los Muchachos on the island of La Palma, Canary 
Islands. One candidate, ULAS J145441.41 + 123556.6 (ID=27), 
was observed on 3 September 2011 (Table |3]l. Weather condi- 
tions were non photometric with high level clouds and seeing 
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Table 2. Revised proper motions and reduced proper motions 
along with the values derived from the VO cross-match 



ID 



/'total 



Hr 



"/yr 



"/yr 



"/yr 



mag "/yr 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



-0.077±0.005 
0.004±0.005 
-0.244±0.001 
-0.022±0.005 
-0.276±0.002 
-0.070±0.002 
-0.147±0.003 
-0.018+0.006 
-0.429±0.006 
-0.141±0.006 
0.040±0.007 

12 -0.026±0.006 

13 0.103+0.007 

14 0.046+0.006 

15 -0.017+0.005 

16 -0.181+0.003 

17 -0.332+0.007 

18 -0.229+0.003 

19 0.063+0.010 

20 -0.280+0.004 

21 -0.193+0.007 

22 -0.525+0.007 
-0.134+0.005 
-0.332+0.006 
-0.352+0.010 
-0.052+0.003 
-0.307+0.010 

28 -0.016+0.008 

29 -0.007+0.006 

30 — + — 

31 — + — 

32 -0.012+0.004 



23 
24 
25 
26 
27 



0.115+0.005 
-0.103+0.005 
0.190+0.002 
0.246+0.004 
0.024+0.001 
0.773+0.001 
0.087+0.004 
0.032+0.006 
0.266+0.005 
0.034+0.006 
0.008+0.008 
-0.029+0.006 
0.069+0.006 
-0.005+0.006 
0.067+0.005 
0.008+0.005 
0.080+0.007 
0.066+0.003 
0.080+0.010 
-0.177+0.003 
-0.259+0.008 
-0.407+0.007 
-0.082+0.004 
-0.120+0.006 
-0.042+0.012 
-0.121+0.004 
-0.104+0.010 
0.015+0.010 
0.008+0.006 



0.139+0.007 
0.103+0.007 
0.309+0.002 
0.247+0.006 
0.277+0.002 
0.776+0.003 
0.171+0.005 
0.037+0.009 
0.504+0.008 
0.145+0.008 
0.041+0.011 
0.039+0.008 
0.124+0.009 
0.046+0.008 
0.069+0.007 
0.181+0.006 
0.341+0.010 
0.239+0.005 
0.101+0.141 
0.331+0.005 
0.323+0.011 
0.664+0.010 
0.157+0.007 
0.353+0.009 
0.355+0.016 
0.132+0.005 
0.324+0.014 
0.022+0.013 
0.011+0.009 



22.215 
22.659 
21.405 
23.050 
21.979 
26.180 
23.662 
21.522 
21.327 
23.127 
20.438 
20.290 
22.681 
21.554 
21.707 
23.249 
23.265 
23.409 
22.919 
22.162 
22.340 
26.610 
21.062 
22.223 
23.256 
22.475 
21.827 
19.136 
18.485 



the mean flat field (also bias subtracted). Second, we extracted a 
one-dimensional spectiTim interactively by choosing the appro- 

priate width for the aperture width. Then, we used the arc lamps 

//(VO) Hr(VO)to calibrate our spectra in wavelength with an accuracy better 
than 0.2 A. The flux calibration of the s pectrum was conducted 



mag 



0.193+0.005 0.193+0.007 22.904 



0.137 
0.135 
0.283 
0.253 
0.270 
0.783 
0.166 
0.297 
0.500 
0.162 
0.614 
1.008 
0.803 
0.273 
0.162 
0.189 
0.344 
0.255 
0.303 
0.290 
0.294 
0.673 
0.178 
0.355 
0.320 
0.188 
0.312 
0.326 
0.733 
1.210 
1.129 
0.202 



22.177 using BDh-174708 dLafliam etal.ll2002h as spectrophotometric 
23.240 standard star. The final optical spectrum, covering the 500-900 
21.21 1 jjjj^ wavelength range and normalised at 750 nm is shown in Fig. 
23.106— 5 6 e 

21.925 
26.201 
23.592 
26.044 



4.2. F0RS2 spectroscopy 



22.514^2 conducted low-resolution (R~350) optical (600-1010 nm) 
23.364 spectroscopy with the visual and near UV FOcal Reducer and 
2j^^^ low dispersion Spectrograph FORS2 (lAppenzeller et al.l Il998h 
26 739 '■^^ ESO VLT Antu unit in Paranal, Chile. We observed 

25 421 22 subdwarf candidates identified in the UKIDSS LAS DR5 
23.556 SDSS DR7 cross-match (Table [3]). All observations were 
23.346 conducted in service mode under programs 084.C-0928A and 
23.284 084.C-0928B. The requested conditions, grey time, thin cirrus 
23.545 acceptable, and seeing less than 1.4 and 1.0 arcsec for the bright 
25.306 and faint targets, respectively, were generally satisfied (TableO. 
21.877 objects were observed at parallactic angle. 

22 133 r o 

26 638 FORS2 instrument is equipped with two 2048x4096 

21 334 CCDs with pixels of 15u m working in the 330-1100 nm 

22 236 r^nge ('A ppenzeller et al.ll99 8V We employed the grism 150-1-27 
23.030 with the order blocking filter OG590 with the standard resolution 
23.248 of 2.07 nm per pixel and a slit of 1 .0 arcsec to achieve a spectral 
21.742 resolution of ~175 at 720 nm due to the 2x2 binning. The expo- 
24.992 sure time was scaled according to the brightness of the target in 
27.602 the Sloan r and / filter sin order to achieve a minimum signal-to- 
2^-1^2 noise of 20. The faintest sources, observed several times along 

the slit, have lower signal-to-noise because they are faint even for 



27.983 
23.005 



around 0.7-0.9 arcsec. The object was observed around UT = 
21h at parallactic angle under high airmass, starting at 1.5 and 
finishing around 1.9 (Table O. 

The ALFOSC spectrograph is equipped with a 2048 x 2052 
pixel back-illuminated CCD42-40 charge coupled device. We 
employed the grism number 5. The total exposure time for 
ULAS J145441.41-H123556.6 was divided into three on-source 
integrations of 800 seconds with a slit width of 1 arcsec to 
achieve a resolution of 1 .55 nm at 700 nm. An internal flat-field 
was obtained immediately after the target to remove as much 
as possible the fringing at long wavelength, the effect being of 
the order of 18% at 800 nm. Unfortunately fringing is clearly 
visible beyond 800 nm because the internal flat field was ob- 
tained after the last of the three individual exposures taken for 
the target. Bias frames were observed in the afternoon before the 
beginning of the night. He, Ne, and Ar arc lamps were also ob- 
tained immediately after each exposure to calibrate our target in 
wavelength. 

The data reduction of the NOT/ ALFOSC spectra was en- 
tirely carried out under the IRAfQ environment. First, we sub- 
tracted the average of all bias frames from the raw science ex- 
posures and then divided by the normaUsed response function of 



* IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation 



a 8-m class telescope. Dome flat fields, bias frames, arc lamps, 
and spectrophotometric standard stars were observed every night 
as part of the ESO calibration plan. 

As for NOT/ALFOSC, the data reduction of the VLT FORS2 
spectra was entirely carried out under the IRAF environment. 
First, we cut the 2D images to select the interesting part of the 
spectrum, from ~600 to 1000 nm roughly. Then, we subtracted 
the average of all bias frames from the raw science exposures 
and from the median-combined dome flat. Afterwards, we di- 
vided the science frame by the normalised flat field using the 
mean value over the entire dome flat frame. Later, we extracted 
a one-dimensional spectrum interactively by choosing the ap- 
propriate value for the aperture width. Then, we used the arc 
lamps with Helium, HgCd, and Argon to calibrate our spectra 
in wavelength with an accuracy of the order of 0.4-0.6 A rms. 
Finally we calibrated the ID spectra with spectrophotometric 
standard stars observed on the same night as the target. Spectra 
have been normalised at 750 nm and are not coiTected for the 
telluric band around 760 nm. The VLT FORS2 spectra of the 1 1 
candidates classified as subdwarfs are displayed in Fig. |4] The 
sources classified as extreme subdwarfs (7 objects) and ultra- 
subdwarf (2 objects) are shown in Fig.|5]while the seven sources 
classified as solar-metallicity M dwarfs are displayed in Fig.|6] 
The two known subdwarfs previously reported in the literature 
are included in these figures. 

4.3. SDSS spectroscopy 

The SDSS spectroscopic database represents an invaluable 
source of good quality spectra for a large number of astronomi- 
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Table 3. Log of the spectroscopic observations 



R.A. 


Dec 


Instr 


Date 


# 


ExpT 


Airm 


Seeing 


hh:mm:ss.ss 


Q .1 .ft 




yyyy-mm-dd 




sec 




" 


00:45:39.97 


+ 13:50:32.7 


FORS 


2009-11-10 


3 


1170 


1.286 


0.74 


01:28:30.89 


+ 13:45:07.4 


FORS 


2009-11-04/15 


4 


2640 


1.286 


0.62 


01:50:34.33 


+ 14:20:02.4 


FORS 


2009-10-30 


3 


1170 


1.494 


0.59 


03:33:50.84 


+00:14:06.1 


FORS 


2009-10-30 


3 


870 


1.546 


0.80 


08:41:53.89 


+02:06:15.1 


FORS 


2009-12-21/25 


4 


2640 


1.121 


0.81 


08:55:00.12 


+00:02:04.1 


FORS 


2009-12-25 


1 


660 


1.122 


0.65 


10:01:26.29 


-01:34:26.6 


FORS 


2010-01-08 


2 


1330 


1.092 


0.86 


10:07:43.96 


-02:28:30.0 


FORS 


2009-12-17/24 


4 


2644 


1.150 


0.83 


10:16:13.89 


+01:13:11.4 


FORS 


2009-12-24 


2 


1320 


1.161 


0.56 


11:58:26.62 


+04:47:46.8 


FORS 


2009-12-21 


1 


1320 


1.468 


0.73 


12:02:14.62 


+07:31:13.8 


FORS 


2010-01-23 


1 


390 


1.222 


0.77 


12:15:08.37 


+04:02:00.5 


FORS 


2010-02-12 


3 


1170 


1.183 


0.87 


12:21:45.28 


+08:04:04.4 


FORS 


2010-01-08/09 


4 


2644 


1.280 


0.60 


12:42:34.62 


+ 14:33:06.2 


FORS 


2010-03-07 


1 


390 


1.306 


0.53 


12:44:25.90 


+ 10:24:41.9 


FORS 


2010-01-23 


1 


390 


1.237 


0.95 


12:46:21.90 


+04:43:09.9 


FORS 


2010-03-07 


1 


390 


1.181 


0.57 


13:17:05.66 


+09:10:16.9 


FORS 


2010-03-07 


1 


390 


1.239 


0.52 


14:18:06.71 


+00:00:35.5 


FORS 


2010-02-13 


3 


1170 


1.119 


0.71 


15:12:11.64 


+06:42:51.3 


FORS 


2010-03-07 


4 


2640 


1.510 


0.57 


15:43:31.93 


+02:45:37.8 


FORS 


2010-03-07 


8 


5280 


1.333 


0.63 


23:33:59.39 


+00:49:35.2 


FORS 


2009-11-10 


3 


1170 


1.130 


1.14 



14 



12 



10 




14:54:41.42 +12:35:56.7 NOT 2011-09-03 3 2400 1.700 0.80 



01:33:46.25 
02:05:33.75 
08:43:58.50 
12:02:14.62 
12:36:59.43 
12:56:35.91 
13:17:05.66 



+ 13:28:22.4 
+ 12:38:24.1 
+06:00:38.6 
+07:31:13.8 
-00:21:58.2 
-00:19:44.9 
+09:10:16.9 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



2000-12-01 

2000- 12-01 

2003- 03-10 

2004- 03-25 

2001- 02-01 
2001-03-26 
2006-04-25 



5400 
4500 
2400 
3000 
3600 
2702 
4900 



1.143 
1.143 
1.161 
1.133 
1.209 
1.214 
1.216 



2.82 
2.82 
1.81 
1.41 
2.87 
2.28 

2.0^ ig.4, 




600 



700 



800 

Wavelength (nm) 



900 



1000 



cal sources. These spectra, covering the 380-940 nm wavelength 
range with a resolution of ~2000 are publicly available from 
the Sloan webpage. They are wavelength and flux calibrated and 
also corrected for telluric absorption bands. The total exposure 
times consist of several integrations of 900 seconds. We found 
seven of our 32 ultracool subdwarf candidates in this database, 
two of them (ULAS J131705.66+091016.9 (ID=25) and ULAS 
J120214.62+073113.8 (ID=17)) re-observed with VLT FORS2 
because of their poorer quality (Table |3). The other five SDSS 
spectra have good signal-to-noise and we analysed them in a 
similar manner as the VLT FORS2 spectra (see Sect. |5}. The ob- 
serving date, number of exposures, airmass and best 80% seeing 
are Usted at the bottom of Table |3] 



5. New ultracool subdwarfs 

In this section, we assign spectral types to the new subdwarfs, 
measure their radial velocities, estimate their spectroscopic dis- 
tances, discuss their mid-infrared properties, search for wide 
companions, and present a preliminary estimate of their surface 
density. 



5A. Spectral types 

We employed two independent but complementary methods to 
assign spectral types to our new ultracool subdwarfs. 

The adopted classication for M-type subdwarfs and extreme 
subdwarfs relies on the scheme proposed by Gizis (1997) and 
extends to spectral types sdM7 and esdM5.5, respectively. This 
scheme is based on the strength of CaH (temperature index) 



Low -resolution optical (600-1000 nm) spectra of con- 
firmed subdwarfs obtained with VLT FORS2. Spectra are or- 
dered by increasing spectral type and offset by a constant value 
for clarity. Overplotted in red are known template subdwarfs 
from the SDSS spectroscopic database (600-940 nm) except for 
the top 3 templates which come from the IRTF/SpeX library. 



and TiO (metallicity index) bands. An extension to later spec- 
tral types (u p to sdM9) has been more recently proposed by 
iLepine et al] ( f2003a) . 

We measured the fou r spect r al indi ces (Ti05, CaHl, CaH2, 
and CaH3) defined by 'Gizis' (1997^ and later updated by 
Lepine et al. (2007) to distinguish ultra-subdwarfs, extreme sub- 
dwarfs, and subdwarfs from solar metallicity M dwarfs based on 
the strength of the CaH and TiO absorption bands (Fig. |7]i. The 
spectral indices and their associated spectral types (quoted to the 
nearest decimal) derived for each of the confirmed subdwarfs are 
fisted in Tableland plotted in Fig.EHE] 

The SDSS spectroscopic database provides optical spectra 
over the 600-940 nm wavelength range for a large number of 
metal-poor dwarfs with spectral types based on two indepen- 
dent classification schemes: the Hammer scheme discussed by 
ICovev et"a n ('2007!) an d the u pdated subdwarf classification pro- 
posed by Lepine et al.l (l2007h . The former distinguishes between 
solar-metallicity M dwarfs and subdwarfs whereas the latter pro- 
vides an accurate classification for metal-poor dwarfs with sub- 
classes. Therefore, we downloaded from the Sloan spectroscopic 
database the brightest object of each subclass between and 9 
for subdwarfs, extreme subdwarfs, and ultra-subdwarfs. These 
"templates" were used to assign visually spectral types to our tar- 
gets with an uncertainty of half a subclass (or better). We should 
mention that the SDSS templates have been corrected for telluric 
absorption whereas our spectra were not, resulting in a differ- 
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600 700 800 900 1000 

Wavelength (nm) 

Fig. 5. Low -resolution optical (600-1000 nm) spectra of con- 
firmed extreme subdwarfs and ultra-subdwaifs obtained with 
VLT FORS2. The second spectiTim from bottom is from 
NOT/ALFOSC. Spectra are ordered by increasing spectral type 
and offset by a constant value for clarity. Overplotted in red are 
known template extreme subdwarfs and ultra-subdwarfs from 
the SDSS spectroscopic database (600-940 nm). 




600 700 800 900 1000 

Wavelength (nm) 



Fig. 6. Low -resolution optical (600-1000 nm) spectra of pho- 
tometric candidates classified as solar-metallicity M dwarfs. 
Spectra are ordered by increasing spectral type and offset by a 
constant value for clarity. Oveiplotted in red are known M dwaif 
templates downloaded from the SDSS spectroscopic database. 
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Fig. 9. Comparison of the spectral types obtained from the vi- 
sual comparison with metal-poor "templates" and spectral t ypes 
derived from the spectral indices defined by IGizisI (1 19971) and 
Lepi ne et al. (2007 ). Subdwarfs, extreme subdwarfs, and ultra- 
subdwarfs are marked as squares, circles, and triangles, respec- 
tively. Metal-poor dwarfs with a discrepant spectral types de- 
rived from both classification methods (comparison with tem- 
plates vs spectral indices) are displayed with asterisks (e.g. sdM 
vs esdM or usdM vs esdM, etc). Typical uncertainties on the 
spectral types derived from both methods are half a subclass. 

ence between them around the position of the O2 telluric line 
around 760 nm. We note that the sdM9.5 template was taken 
from the SpeX library^ because no sdM9.5 template was found 
in the SDSS spectroscopic database. The resulting fits (red) to 
our spectra (black) are presented in Figures H]to |6] 

The comparison between the spectral types inferred from 
spectral indices and from the direct comparison with "templates" 
is shown in Fig. |9] We find that the spectral types derived from 
spectral indices tend to under-estimate the spectral type (over- 
estimate the effective temperature). Therefore, we adopted the 
direct and visual comparison with "templates" to assign spec- 
tral types to our targets because it provides a more accurate 
classification. We note that the spectral indices are not so reli- 
able to classify subdwarfs bec ause they rely on a n arrow wave- 
length range as pointed out by iLepine et al.l (l2007h and also de- 
pend strongly on the resolution of the spectra (Table|4). The dif- 
ferences between both classifications are listed in Table |4] and 
shown in Fig. |9] where asterisks represent confirmed metal-poor 
dwarfs with discrepant classes (not only spectral types) derived 
from the spectral indices and the direct comparison with tem- 
plates, suggesting some shortcomings in the spectral types de- 
rived solely on indices. Our sample of 20 low-mass stars with 
low-metal content consists of nine subdwarfs, seven extreme 
subdwarfs, two ultra-subdwarfs, and two L subdwarfs. 

The remaining seven sources in our spectroscopic sample are 
M dwarfs with spectral types ranging from M4 to M7. We be- 
lieve that the contamination of our sample by these M dwarfs 

^ http://web.mit.edu/ajb/www/browndwarfs/spexprism/ 
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Table 4. Coordina tes (in J2000), spectral indices, spectral types determined following the definitions by iGizisI (Il997h and 
iLepine et alJ (l2007 l) for the new subdwarfs. The last column gives the adopted spectral types derived from direct comparison with 
spectral templates. If a target appears twice, the first line corresponds to the FORS2 spectrum while the second is the SDSS spectrum. 
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come from large error bars on their Sloan positions, leading to 
spurious proper motions placing them at the bottom of the re- 
duced proper motion diagram mimicking subdwarf candidates. 
We remark that a revision of proper motions as that outlined in 
Section [3] likely reduces the contamination level by a factor of 
2-3. None of these M dwarfs exhibit Ha in emission at the reso- 
lution of our spectra (R~150), suggesting that they are older than 
typical M dwarfs in the solar vicinity. West et al.l (120081) reported 
Ha in emission for 50% of the sources or more at spectral types 
later than M4, with equivalent widths larger than 1 A. These au- 
thors also discussed the decrease in activity with age and scale 
height, suggesting that the M dwarfs contaminating our sample 
may be either older than the average low-mass stars in the solar 
neighbourhood and may be located at a higher scale height. The 
lack of Ha in emission also places a lower limit on the ages of 
these M dwarfs, 5 and 8 Gyr for M5 and M7 dwarfs, respectively. 

5.2. Radial velocities 

In this section we tentatively compute the radial velocities of our 
targets in spite of the low-resolution of our spectra. We used two 
different methods. 

First, we computed the offsets between lines resolved in our 
FORS2 spectra and the centroids of several atomic lines whose 
accurate positions can be found on the webpage of the National 
Institute of Standards and TechnologyQ: Cal at 6572. 18 A, Til 
at 8434.94 A, Nal doublet at 8183.25 and 8194.79 A, Csl at 
8542 09 and 8943.47 A, an d Call at 8542.09 A (see also Table 
2 of iBurgasser et al.l l2009l) . The uncertainties on these offsets 
are typically of the order of 1 A, corresponding to typical un- 
certainties of 35-50 kms"'. We focused only on the Cal, Nal 
doublet, and Call lines for the NOT spectrum and the SDSS 
spectra. We derived the observed radial velocities by multiplying 

* |http://physics.nist.gov/asd3 1 



those offsets by the speed of light and dividing by the wavelength 
{c*AA/A). Observed velocities were converted into heliocentric 
velocities by computing the Earth's rotation, the motion of the 
Earth's center about the Earth-Moon barycenter, and the motion 
of the Earth-Moon barycenter about the center of the Sun using 
the RVCORRECT routine in IRAF In the fifth column of Table [f] 
we provide our heliocentric radial velocities measured with this 
method. 



Independently, we computed radial velocities by cross- 
correlating our optical spectra with the optical spec trum of the 
sdL3.5 subdwarf SDSS J125637. 13-022452.4 ( Bur gasser et akl 
2009) using the task FXCOR under the IRAF environment. This 
subdwarf has a well-measured radial velocity of -130 km/s and 
we decided to use it as a template despite having a spectral type 
later than our confirmed s ubdwarfs. We transform ed the vacuum 
wavelength published by iBurgasser et al.l (l2009l) into air wave- 
length and degraded its resolution to match our observations. 
We cross-matched the full spectrum of our sdM9.5 and sdL sub- 
dwarfs with Burgasser's template because of the similarities in 
the shape of the spectra and spectral types but focused only on 
the Nal doublet region for the earlier subdwarfs in our sample. 
We cross-matched the spectra using two different functions: a 
gaussian fit to the cross-correlation function peak or simply the 
central value of the peak. Both method led to differences smaller 
than the uncertainties output by FXCOR. The radial velocities in- 
ferred by this method are listed in Table |5]and have typical un- 
certainties of 40-50 km s"' measured from the dispersion of the 
measurements from different part of the spectra. Radial veloci- 
ties derived from both methods usually agree within the uncer- 
tainties, except for ULAS Jl 15826.62+044746.8 (ID=16) and 
ULAS J12:36:59.43-00:21:58.2. In this case, we favour the ra- 
dial velocities derived from the direct comparison with SDSS 
J125637. 13-022452.4. 
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Fig. 7. Top left: CaHl vs Ti05 indices for our new ultracool subdwarfs. Top right: CaH2 vs Ti05 indices. Bottom left: CaH3 
vs Ti05 indices. Bottom right: Sum of CaH2 and CaH3 vs Ti05 indices. Subdwarfs, extreme subdwarfs, ultra-subdwarfs, old 
solar-metallicity M dwarfs from our sample are marked as filled squares, circles, triangles, and diamonds, respectively. Numbers 
d enomin ating our new discoveries follow the order by right ascension from Table|4] Open symbols are known subdwarfs discussed 
m lGizislfl997h . The small coloured dots in the bottom right plot represent sources with SDSS spectroscopy classified as subdwarfs 
(red), extreme subdwarfs (green), and ultra-subdwarfs (blue). For the subdwarfs with two spectra from FORS2 and SDSS, we plotted 
only the indices derived from the VLT FORS2 spectra. Typical uncertainties on the spect ral indices a re of the order of 0.1 (cross at 
the bo ttom of each plot). These plots follow the standard figures presented in Figure 1 o f|G izisi (1 19971) and Figure 3 of iLepine et al.l 
(l2007h . 



5.3. Spectroscopic distances 

In this section, we estimate the spectroscopic distances by com- 
paring our new discoveries with subdwarfs of similar spectral 
types with known trigonometric parallaxes. 

We looked for subdwarfs with parallaxes whose spec- 
tral types are sdM5, sdM6, sdM6.5, sdM7, sdM7.5, and 



sdM9.5. We considered the following sources as subdwarfs 
with known distances to derive spectroscopic distances for 
our new ultracool subdwar fs: LP 807-23 (sdM5.0 ; J = 12.92 
mag; d = 28.17-31.74 pc: Ivan Altenaetan IT991. LHS 1074 
(sdM6; J ^ 14 68 ma g; d = 85.7+17.1 pc: ISaUm & Gouldl 
l200l iRiaz et alJl2008h. LHS 1166 rsdM6.5: J = 14.26 mag; 

in [199: 



73-89 pc; Ivan Altena et all 119951) . LP 440-52 (sdM7; 
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Fig. 8. Spectral indices as a function of the adopted spectral types for our new ultracool subdwarfs. Typical uncertainties on the 
spectral indices are of the order of 0.1 (cross at the bottom of each plot). Uncertainties on the spectral types are 0.5 subclass. Top 
left: Adopted spectral types vs Ti05 index. Top right: Adopted spectral types vs CaHl index. Bottom left: Adopted spectral types 
vs CaH2 index. Bottom right: Adopted spectral types vs CaH3 index. 



/ = 13.19 mag; d = 34.4-36.1 pc; IvanAltena et"al] Il995h . 
LSRJ20362 1.86+ 505950. 3 (sdM7.5; / ^ 13.628 mag; d 
= 43.7-49.2 pc; iLenine et al.l l2002t ISchilbach et all l2009h . 
LSR J 142504. 8 1 +7 102 10.4 (sdM8; J ^ 14.828 mag ; d = 
75.4-89.9 pc; iLepine et al.l 120031; iBurgasser et all 120081: 
2009h . and SSSPM101307. 34-135620.4 
mag; d = 45.0-54.6 pc; IScholz et all 



ISchilbach et al.1 
(sdM9 .5; J = 14.637 



I2004al; ISchilbach et"!!] l2009l) . After applying the standard 
transformations using the 7-band magnitudes of our new 
discoveries, we assigned mean distances between 88 and 628 
pc (see Table |6]l, assuming that they are single. Typical eiTor 
bars on the spectroscopic distances are 20-25% taking into 



account the uncertainties on the trigonometric parallaxes of the 
templates. For the two L subdwarfs, we used the trigonometric 
parallax of the sdM9.5 template to place upper limits on their 
distances. Both objects very likely lie within 100 pc unless they 
are binaries (Table|6]l. 

In addition to the aforementioned subdwarfs, we were able 
to assign a spectroscopic distance for our esdM5 extreme sub- 
dwarf using LHS515 as template (e sdM5; J - 13.64 mag; d = 
42.6-64.5 pc van Altena et al ] IT995h . Hence, we derive a spec- 
troscopic distance of 257 pc with a probable range of 214-323 
pc for ULAS J124621.90+044309.9 (ID=23; Table |6]l. For the 
other extreme subdwarfs we are unable to assign spectroscopic 
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Table 5. Coordinates (in J2000), spectral types, and radial veloc- 
ities (expressed in km s"') derived in two different ways: first by 
measuring the offsets of well-known lines comparing with their 
centers (typical uncertainties of 35-50 kms"') and, second, by 
cross-coiTelating the Nal doublet region of our optical spectra 
using SDSS J125637. 13-022452.4 (sdL3.5) as template. 
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-163.9 


-192.9 



Table 6. Coordinates (in J2000), spectral types, and spectro- 
scopic distance estimates (in pc) for our new subdwarfs. 



ID 


R.A. 


Dec 


SpT 


Distance 




hh:mm:ss.ss 






pc 


1 


00:45:39.97 


4-13:50:32.7 


sdM7.5 


354 (334-376) 


3 


01:33:46.25 


-^13:28:22.4 


esdM6.5 


124° 


4 


01:50:34.33 


4-14:20:02.4 


sdM6.0 


303 (243-364) 


5 


02:05:33.75 


4-12:38:24.1 


sdMS.O 


132 (121-146) 


6 


03:33:50.84 


4-00:14:06.1 


sdLO.O 


<97 (88-107) 


9 


08:43:58.50 


4-06:00:38.6 


usdM5.5 


no template 


10 


08:55:00.12 


4-00:02:04.1 


sdM6.0 


524 (419-628) 


13 


10:01:26.29 


-01:34:26.6 


sdM5.0 


464 (438-494) 


16 


11:58:26.62 


-(-04:47:46.8 


sdM9.5 


116(105-128) 


17 


12:02:14.62 


4-07:31:13.8 


esdM7.0 


237° 


18 


12:15:08.37 


4-04:02:00.5 


sdM7.0 


281 (274-288) 


20 


12:36:59.43 


-00:21:58.2 


esdM5.5 


171° 


21 


12:42:34.62 


4-14:33:06.2 


usdM5.0 


no template 


22 


12:44:25.90 


4-10:24:41.9 


sdLO.5 


<104 (95-115) 


23 


12:46:21.90 


4-04:43:09.9 


esdM5.0 


257 (214-323) 


24 


12:56:35.91 


-00:19:44.9 


esdM6.0 


132° 


25 


13:17:05.66 


-h09:10:16.9 


esdM7.0 


214° 


25 


13:17:05.66 


-h09:10:16.9 


esdM6.5 


242° 


26 


14:18:06.71 


4-00:00:35.5 


sdM6.5 


464(424-517) 


27 


14:54:41.42 


4-12:35:56.7 


esdM5.5 


150° 


32 


23:33:59.39 


4-00:49:35.2 


sdM6.0 


399 (319-478) 



° the range in distances is not listed because the templates used do not 
have trigonometric parallaxes. 



distances because no object with similar spectral types has 
known trigonometric parallax in the literature. Instead, we used 
LHS 2096 (esdM5.5; J = 13.99; d = 56.10 pc; 'Lepine & Shara' 
2001 . LHS2023 (esdM6; J = 14.91; d = 73.9 pc; Riaz et al. 
20081) . LSRJ0822+1700 (esdM6.5; / = 15.517; d = 106 pc; 



Lepine et al.ll2003d). and A PMPM0559 (esdM7; J = 14.887; d 
— 70 pc; Schweitzer et al.i ri999) to infer tentative (mean) dis- 
tances. We do not quote uncertainties for these sources because 
of the (akeady) very uncertain distances of the templates used. 
We do not provide a distance for the two ultra-subdwarfs in our 
sample because none has trigonometric parallaxes published in 
the literature. 

5.4. Mid-infrared pliotometry 

We cross-matched our list of new ultracool subdwarfs with the 
Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) 
data release which took place in April 201 1. We found two sub- 
dwarfs with mid-infrared photometry using a matching radius of 
6.5 arcsec, the spatial resolution of WISE. However, we ensured 
that no other source was detected in WISE because of the large 
proper motions of subdwarfs. The WISE photometry of these 
two subdwarfs at 3.4 (Wl) and 4.6 (W2) microns is reported in 
Table |7] and plotted in Fig.[TOl these two sources are undetected 
at longer wavelengths. 

We compared the infrared colours of our new u ltracool sub- 
dwarfs to the sample of Kirkp atrick et al] (1201 Ih drawn from 
the DwarfArchives.org webpageQ in order to identify any trend 
that may help future searches for metal-poor brown dwarfs. In 
the case of ULAS J015034.33-H142002.4 (ID=4; sdM6), we 
find colours similar to normal M6 dwarfs in H-W2 - 1.07, 
7-W2=0.76, and K-Wl although on the red side of the distri- 



bution. Other colours, including /-Wl = 0.76, W1-W2 = 0.79, 
K-Wl - 0.92 mag differ to those of normal M6 dwarfs which 
span the following ranges 1.0-1.4, 0.2-0.4, and 0.3-0.6 mag, 
resulting in a deviation of ~lcr (Fig. [TOl Table |7]). The other 
subdwarf in our sample, ULAS145441.42 4-123556.7 (ID=27; 
sdM5.5), exhibits similar colours to M dwarfs in all combina- 
tions of colours at odds with ULAS J015034.334-142002.4 (Fig. 
[Tot . We repeated the process with a larger number of known 
subdwarfs and did not spot any obvious trend with decreasing 
metallicity, suggesting that ULAS J015034.334- 142002.4 may 
be peculiar, or a multiple source, or more likely the WISE pho- 
tometry has underestimated uncertainties. 

We also checked the Spitzer archive to extract additional in- 
formation on the mid-infrared properties of some of our can- 
didates but none of them was found from the Spitzer public 
database. 



5.5. Search for wide companions 

In this section we looked for wide companions brighter than each 
of our subdwarfs within a radius of 10 arcmin. The idea is to find 
potential primaries with distances, metallic ities, and (possibly 
ages (often refered to as benchmark objects; iPinfield et al.ll200 
determined with a higher precision than our new subdwarfs. 
We used larg e-scale surveys wi th accurate prop er motions: the 



USNO -Bl (Monet et alJ l2003h. the UCAC3 (IZacharias et all 



http://spider.ipac.caltech.edu/staff/davy/ARCHIVE/index.shtml 



i2010h . and the PPMXL dRoeser et alJl2010h catalogues. We se- 
lected only bright sources with /usno ^18 mag and proper mo- 
tions in right ascension and declination within 30% of the mea- 
sured motion of our targets derived from the UKIDSS LAS DR5 
vs SDSS DR7 cross-match. 

We found one potential wide companion to ULAS 
J233359.394-004935.2 (ID=32) in USNO-Bl and PPMXL at 
a distance of about 5 arcmin on the sky. This potential com- 
panion, USNO J233350.274-005342.9, has a proper motion in 
right ascension and declination of (149, -75) and (154, -78) 
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Table 7. Coordinates (in J2000), spectral type, 7-band magnitude, and mid-infrared photometry with the associated error bars (3.4 
and 4.6 microns) for the two ultracool subdwarfs covered by the WISE mission. The signal-to-noise ratio of the photometry quoted 
by the WISE catalogue is indicated in brackets. The photometry in the WISE 12 and 22 micron bands is not included because the 
signal-to-noise ratios are less than 3. 



ID 


R.A. 


Dec 


SpType 


J 


3.4yum (SNR) 


4.6yum (SNR) 




hh:mm:ss.ss 






mag 


mag±eir () 


mag+eiT () 


4 


01:50:34.33 


+ 14:20:02.4 


sdM6.0 


17.424 


16.667±0.131 (8.3) 


15.877+0.230 (4.7) 


25 


14:54:41.42 


+ 12:35:56.7 


sdM5.5 


16.121 


15.136+0.044 (24.4) 


14.837+0.090 (12.1) 



1.0 

0.8 h 

(M 0.6 

I 

^ 0.4 

0.2 



0.0 



2 4 6 

M Spectral Type 



10 



This source 



Fig. 10. WISE colours as a function of spectral type for 

known M dwarfs (open squares) from the sample published by 

iKirkpatricket all (12011 '). Our two subdwarf candidates are high- iLepine & 
lighted with large filled circles. 



5.6. Notes on individual objects 

In this section we give additional details on a few specific candi- 
dates identified in our search for ultracool subdwarfs. 



5.6.1 . Spectra in the SDSS DR7 spectroscopic database 

Some photometric candidates in our sample have spectra pub- 
licly available in the SDSS DRV spectroscopic database (Table 
O. As noted in the caption of Table [1] in Appendix, all of them 
were included in the sample of West et al. (2008) a nd classified 
as ear ly-M dwarfs using the Hammer classification (ICovev et al.l 
|2007 |). However, the Sloan spectra clearly look like subdwarf 
with strong CaH and TiO bands typical of low-metallicity 
M dwai-fs (Figs. HHU. We should mention that LP 468-277 
(01:33:46.25-Hl3: 28:22.4) was iri c luded in the catalogue of 
Northern stars of iLepine & Sharal (l2005h but no spectral type 
was derived. Finally, we recovered SDSS J020533.75H-123824. 1 
(ID=5) classified as sdM7.5 by Lepine & Scholz (.2008,) and re- 
classified as sdM8 in this paper based on the spectral type pro- 
vided by the SDSS DR7 spectroscopic database. 



5.6.2. ULASJ1 45441 .45+1 23557.6 (ID = 27) 



ce IS part o 
Sharal (120051) 



art o f the catalogue of Northern stars of 
but no spectral type was derived. Our 



proper motion derived from the LAS and SDSS DR7 epochs 
(0.31 arcsec/yr) is in good agreement \ yith th e value of 0.321 
arcsec/yr published bv I Lepine & Sharal (l2005l) . We do not have 
spectroscopic follow-up for this source yet. 



mas/yr in PPXML and USNO-Bl, respectively, in agreement 
with the values of (136, -71) mas/yr quoted by Sloan. The dif- 
ferences in the proper motions in light ascension and in declina- 
tion between our subdwarf and its potential companion are about 
23-27% and 10-15%, respectively. No spectrum is available in 
the SDSS spectroscopic database. The SDSS/ of this object is 
19.118+0.021 mag, roughly 1 mag brighter than our spectro- 
scopic subdwarf. Its optical colours (g - r- 1.785, r-i- 1.429, 
r-z = 2.158 mag) and reduced proper motion (H^ - 21.64 mag) 
satisfy our original selection criteria and suggests that this po- 
tential wide companion may also be metal-poor. However, its 
J - K infrared are redder than our original cut of 0.7 mag, im- 
plying that this object felt out of our sample. At the spectroscopic 
distance of ULAS J233359.39+004935.2 (ID=32) estimated to 
400 pc (Table |6]l, the projected physical separation of the sys- 
tem would be very larger, of the order of 120,000 au. Hence, we 
cannot claim that both objects are gravitationally bound but they 
might have formed in the same cluster or might belong to the 
same moving group. 



5.6.3. ULAS J233359.39+004935.2 (ID = 32) 

This candidate, observed with VLT FORS2, is classified as 
a sdM6 subdwarf. Its proper motion in right ascension and 
declination is (93.9, 23.0) and (81.9, 16.6) mas/yr from the 
2MASS/UKIDSS and SDSS/UKIDSS cross-match. We discuss 
the presence of a possible wide companion five arcmin away in 
Sectionl53] 

This subdwarf is located at ~42 arcsec fr om 2MASS 
J233358.40+005011.9, a LO dwarf reported by IZhang et al.l 
(120101) with a proper moti on of (139.7, 29.5 ) mas/yr reported 
by the PPMXL catalogue (iRoeser et alJ l201Cf). The proper mo- 
tion values in right ascension differ by 50-70%, suggesting that 
both objects are not physically associated. Nonetheless, we in- 
vestigated further the SDSS DR7 spectrum of this possible wide 
companion and compared it to known M9 and M9.5 dwarfs, 
young L dwarf templates (ICruz et al ] 12009"), a sdM9.5 subd- 
warf (Scholz etal. 20 04al) . and our two possible L subdwarfs 
(Fig. [TTT i. The best fit is obtained for a field M9 dwarf of so- 
lar metallicity, LHS 2065 (iKirkpatrick et al.lll99ll) downloaded 
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4 

2 


600 



2MASSJ23335840p00501 1 9 

LHS2065 (dM9.0) 

SSSPM J101 3-1356 (sdM9.5) 




650 



700 750 800 
Wavelength (nm) 



850 



900 



Fig. 11. Comparison of 2MASS J233358.40+00501 1.9 classi- 
fied as LO by Zh ang et all ( 1201 Oh with a known sol a r-meta llicity, 
LHS 2065, classified as M9 by 'Kirk patrick etaP (Il99ll) and a 
sdM9.5 subdwarf dScholz et al...2004al) . 



from Kelle Cruz's webpage (Fig . [TTlPI, earlier by o ne subclass 
compared to the classification bv lZhang et al.l (l2010t) . 

Overall conclusion: these two sources separated by about 42 
arcsec are unlikely to be physically associated. Note that we also 
discussed in the previous section (Sect. l53] l that this object may 
be part of a wide binary system or a old moving group. 

5.6.4. ULASJ033350.84-I-001 406.1 (ID = 6) 

This blight source (J - 16.1 mag) has a large proper motion of 
0.78 arcsec/yr. Its VLT FORS2 spectrum appears slightly red- 
der than ULAS Jl 15826.62 -1-044746.8 (1D=16) which is well 
reproduced by the IRTF/SpeX spectrum of SSSPM 1013-1356 
(sdM9.5; Scholz et al. 2004a; Burgasser 2004). Hence, we pro- 
pose this object as a new sdLO template for future searches be- 
cause it is the first of its class known to date. 

5.6.5. ULAS J1 24425.90-1-1 02441 .9 (ID = 22) 

This bright (J = 16.26 mag) candidate is the coolest object in our 
sample, cooler than ULAS J033350.84-f001406.1 (1D=6; see 
previous section). Hence we propose it as a sdLO.5 template for 
future works on subdwarfs because it is the first of its class at the 
time of writing. Its proper motion also exceeds the 0.5 arcsec/yr 
threshold of high proper motion sources, with a total motion of 
0.67 arcsec/yr. We note that the two new targets classified as L 
subdwarfs occupy the lower right of the reduced proper motion 
diagram (Fig.|2]i, suggesting a steep reddening in the colours of 
subdwarfs in the M/L transition. 



5. 7. Surface density of subdwarfs 

In this section, we provide a tentative estimate of the surface den- 
sity of ultracool subdwarfs i.e. metal-poor dwarfs with spectral 
types later than M5. We found a total of 18 M5-M9.5 subdwarfs 
(and two additional subdwarfs classified as L0-L0.5) over 1343 
square degrees common to UKIDSS LAS DR5 and SDSS DR7 
and imaged in the YJHK filters. We should potentially add one 
source without spectroscopy yet. We derive a density of 0.015- 
0.0 16 >M5 subdwarf s per square degree. 

IWest et al.l (l201ll) identified 70,841 spectroscopic M0-M9 
dwarfs in the 8200 square degrees covered by the SDSS spec- 
troscopic database. However the SDSS spectroscopic follow-up 



of the imaging survey is incomplet e. For this reason , we co nsid- 
ered the photometric sample from 'Bochanski e t al.l (12010 1. their 
Figure 4) and focus on sources with / brighter than 22 mag, 
r-z>2.5, i-z> 0-2, and r-/>0.3 mag, corresponding to >M5 
dwarfs and later. We counted a total number of 653,625 photo- 
metric > M5 dwarfs in 8000 square degrees surveyed by SDSS 
DR^ implying a density of ~82 > M5 dwarfs per square degree. 
We sent a query with th e aforem entioned criteria to the WFCAM 
Science Archive (Ham blv et al.l l2008) to see how many M5 (and 
later) dwarfs in UKIDSS LAS DR5 and SDSS DR7 we could 
recover in order to match those numbers with our search crite- 
ria defined in Sect. |2] We imposed a detection in YJH but not 
in K and requested good quality point sources in addition to the 
constraints on the optical colours. The query returned 113,393 
sources in 1343 square degrees. I f we hmit th e sample to dwarfs 
with z - J greater than 1.4 (West et ani201 ll) . the query returns 
106,746 sources, implying a number of M5 dwarfs (and later) of 
the order of 79.5-84.4 per square degree which is highly con- 
sistent with the numbers derived from the Sloan sample alone. 
This density is ~5000-5700 times higher than the number of ul- 
tracool subdwarfs found in our photometric and proper motion 
search which is broadly consist ent with the 0.2% contribution 
from metal-poor stars quoted by Digbv et al.l (l2003b and the up- 
per limit derived from the SDSS M dwarf sample (ICovev et al.l 
[20081) . We should mention that according to the evolutionary 
models (fBaraffe et al. 1997, 1998), the masses of >M5 dwarfs 
and subdwarfs are similar at ages of Gyr but lower metallicity 
M dwarfs have high effective temperatures. Finally, we should 
point out that we found two ultra-subdwarf and seven extreme 
subdwarfs for 11 subdwarfs, suggesting a fairly quick decrease 
in the numbers of subdwarfs as a function of metallicity. 



6. Summary 

We have presented the outcome of a dedicated photometric and 
proper motion search aimed at finding new ultracool subdwarfs 
in public databases. We identified 32 ultracool subdwarf can- 
didates, 20 of them being confirmed as metal-poor late-M and 
early-L dwarfs by low-resolution optical spectroscopy. We dis- 
covered two new early-L subdwarfs which we propose as spec- 
troscopic templates for future searches because these are the first 
of their subclass. We measured radial velocities for most of the 
new subdwarfs with the cross-match technique. We estimated 
their spectroscopic distances when templates of similar spectral 
types with trigonometric parallaxes were available. We uncov- 
ered seven old M dwarfs contaminating our sample whose ages 
are estimated to >5-8 Gyr due to the lack of Ha in emission. Of 
the 32 candidates, five do not have optical spectroscopy. Only 
one of these five remains a good subdwarf candidate, the oth- 
ers being rejected. We found a contamination of about 30% by 
solar-metallicity M dwarfs in our photometric and proper mo- 
tion search, mainly due to large errors on the Sloan positions 
leading to spurious proper motions affecting the determination 
of the reduced proper motion. We are able to reduce this level 
of contamination by a factor 2 to 3 after revision of the proper 
motion measurements. We also present mid-infrared data from 
WISE for two subdwarfs as well as a search for bright and wide 
common proper motions which led to an extremely wide pair 
very likely not gravitationally bound. 

We intend to expand our search for subdwarfs with upcom- 
ing data releases from UKIDSS to increase the census of low- 
metallicity dwarfs. Moreover, we plan to update our colour crite- 



http://kellecruz.com/M_standards/ 



Numbers kindly provided by John Bochanski and Andrew West 
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ria to optimize future searches and discover even cooler ultracool 
subdwarfs. The main overall scientific goal of this large project 
is to update and extend the current low-metallicity classification 
into the L dwarf (and later T dwarf) regime. 
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2.643' X 1.94' 



N 



2.643' X 1.941' 




Fig. .1. Finding charts for 16 confirmed subdwarfs. Images are in the Sloan z-band and approximately 2.64x1.94 arcmin aside with 
North up and East left. The ID number of the object is indicated at the top left of the image and the scale at the bottom left. 
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Fig. .2. The other 16 candidates identified in this work: the remaining four subdwarfs confirmed spectroscopically are shown at the 
top, the seven solar-metallicity M dwarfs below followed by the five candidates without optical spectra. Images are in the Sloan 
z-band and approximately 2.64x1.94 arcmin aside with North up and East left. The ID number of the object is indicated at the top 
left of the image and the scale at the bottom left. 
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